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Abstract 

We assume that a self-gravitating thin string can be locally de- 
scribed by what we shall call a smoothed cone. If we impose a specific 
constraint on the model of the string, then its central line obeys the 
Nambu-Goto equations. If no constraint is added, then the worldsheet 
of the central line is a totally geodesic surface. 
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I Introduction 



In the framework of general relativity, two theoretical aspects of cosmic 
strings have been mainly studied. The dynamical one in which the equations 
of motion of the infinitely thin cosmic string are governed by the Nambu- 
Goto action and secondly the self-gravitating one in which the straight cosmic 
string in particular is considered as the source of a gravitational field. In this 
case one finds that the asymptotical metric generated by the straight cosmic 
string is a conical metric and the case of a singular line is obtained in a 
limit process. 

The study of line sources in gravitation can be originated from a paper of 
Israel in which the author concluded that there existed " no simple general 
prescription [...] for obtaining the physical characteristics of an arbitrary 
line source" . However, Vilenkin gave the physical meaning of the conical 
singularity as being a cosmic string. Attention has been recently devoted 
in understanding the dynamics of a conical-type line source of the Einstein 
equations [H, |^, ||, ^. The conclusions summarised in |^ are as follows. A 
two-dimensional timelike worldsheet whose points are conical singularities of 
the 4-geometry cannot be in general the dynamical evolution of a Nambu- 
Goto string of arbitrary initial shape. The conical singularity requires that 
the worldsheet is totally geodesic. This restraints the initial shape of the 
string as well as the evolution of each of its points to be a geodesic of the 
4-geometry. 

In this paper, we take up again the question of the dynamics but for a 
self-gravitating extended string. For a straight cosmic string, it is known 
that the exterior metric may be matched with an interior metric having as 
source an energy-momentum tensor with suitable properties. We consider a 
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string of arbitrary shape but we restrict ourselves to the case of a thin tube 
of matter adopting the assumption that the exterior metric is locally the one 
describing a straight cosmic string. The interior metric is constrained from 
this and the energy-momentum tensor could be derived from the metric. The 
aim of the present work is to find the equations of motion of the central line 
of this thin tube of matter in the limit where its radius tends to zero. 

We emphasize that we do not consider directly a self-gravitating line. 
In our method, the conical points are smoothed out on a scale comparable 
to the radius of the string. The central line of this tube sweeps a timelike 
2-worldsheet whose points are perfectly regular. A local coordinate system 
can then be attached to our spacetime by taking the two parameters of 
the worldsheet as the first two coordinates and the other two as geodesic 
coordinates pointing in a direction orthogonal to the worldsheet. This local 
coordinate system affixed to the worldsheet allows the natural introduction 
of the extrinsic curvature and other geometric parameters of the worldsheet. 

The interior metric of the string is essentially characterised by a function 
/(/) = th{l/e) in which / is the radial coordinate whose origin lies on the 
worldsheet and e is a length typical of the thickness of the string. The 
function h is arbitrary, only submitted to certain conditions expressing the 
smoothness of the spacetime on the worldsheet, as well as the matching 
conditions on the boundary between the interior of the string and the vacuum. 

If no other conditions are imposed, we shall prove that the worldsheet 
tends to a totally geodesic surface, i.e. of vanishing extrinsic curvature, 
when e goes to zero. However wc have found another possibility. If we im- 
pose on the function h a specific supplementary constraint then the extrinsic 
curvature no longer tends to zero. Nevertheless the mean curvature tends to 
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zero when e does and thus the worldsheet tends to be locally extremal which 
is precisely the behavior of the Nambu-Goto string. Since the function h is 
a characteristic of the metric in the interior of the string thus the constraint 
imposed is interpreted as a specification of the matter of the string. 

The paper is organised in the following manner. We recall in Sec. II the 
basic results on self-gravitating straight strings with some thickness and we 
describe the formalism that we shall use. In Sec. Ill, a self-gravitating string 
of arbitrary shape is introduced and the coordinate system adopted to its 
study is introduced. In Sec. IV, we expand in powers of 1/e the geometrical 
quantities appearing in the Einstein equations. In Sec. V, taking the limit 
of e going to zero of the Einstein equations on the boundary between the 
string and the vacuum, we obtain the constraints that the worldsheet swept 
by the central line of the string must satisfy. Finally in Sec. VI, we find the 
Nambu-Goto energy-momentum tensor in the zero limit of e, confirming the 
choice of the definitions adopted in Sec. III. 



II Straight strings as smoothed cones 

In a general relativistic context, it is possible to portray a straight string as 
a thin cylinder of matter [|, |ll], |12|. We give some basic features of a 
self-gravitating straight string since they will be used later on in the general 
case of strings of arbitrary shape. 

In the coordinate system [t, z, I, (p) with I > and < < 27r, the 
energy-momentum tensor has the form 

T} = T: = -a{l) T/ = T; = 0</</o (1) 
where Iq is the given radius of the cylinder. The energy density a is a posi- 
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tive regular function of 3?^. Taking into account the Einstein equations, the 
interior metric can be written as 

dsjj^T = -dt^ + dz^ + dl'^ + f{l)d(j)\ 0</</o (2) 

where the positive function / determines a by the formula, 

1 /" 
" = "8^7' 

G being the Newtonian gravitational constant. In order to ensure a regular 
behavior of metric (|^) at / = 0, one must choose / such that 

/(/) ~ / + asP + a^l^ + ■■■ as / ^ 0. (4) 

We can also require that the function / is increasing. The exterior metric 
which can be matched to the interior metric (0) can be expressed in the form 

dsEXT = -dt^ + dz^ + dt^ + sin^ a [l- l^fdcj)^, I > Iq (5) 

where the constants a and Iq are determined from the following matching 
conditions of the metric and its first derivatives at / = /q 

/(/q) = sinQ;(/o— ^o) and /'(/o) = sin a. (6) 
The linear mass density /i of the straight string is defined by 



/i = / af dl d(f). (7) 

JKlo 

By using (^, and (^, it is easy to see that it is related to the angular 
deficit A of metric (|^) by 

A = 27r(l - sin a) = SnGfx. (8) 
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From (H), we note that the exterior metric (|^) is independent of the details 
of the internal string as well as of its radius. 

At this point, it is enlightening to recall the solution for a straight string 
with a =constant. The interior metric is characterised by 

/(0=esin- and a = (9) 

and the exterior metric is again (^. The geometrical interpretation of the pa- 
rameter e will be given below. The matching conditions (|^) give the relation 
/o/e = 7r/2 — a and thereby 

Tn sin a . . 

e={k-h) . 10 

cos a 

We point out that lo/e depends only on a. 

Returning to the general case this particular solution suggests to impose 
in the generic interior metric (Q) the following form 

/(/) = eh{-) (11) 

where is a smooth function and e is a parameter which takes again the 
value ([l0|) . Now the matching conditions (|^) yield simply 

/i(— )=cosa and /i'(— ) = sina. (12) 

For a given function /i, the quotient /o/e depends only on the constant angle 
a, i.e. on the linear mass density /i by virtue of relation (|]). 

One can give a geometrical interpretation. By performing the change of 
coordinate p = I— Iq in metric (^, we recognise the so called conical metric 
representing a cone of half angle a. It is usefuU to give a representation 
in a Euclidean 3-space of the 2-surface t = constant and z = constant for 
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the interior and exterior metrics (cf. fig. 1). It can be visuahsed as a cone 
whose top is cut out at a distance po = of the vortex and replaced 

by an axisymmetric cap which joins the cone tangently along the circle of 
junction of radius po sin a. The coordinate I represents the length of the 
radial geodesic of the 2-surface originated at the top of the cap. It takes the 
value /q at the boundary between the smooth cap and the cone and then it 
continues along a generatrix of the cone. The coordinate is the azimutal 
angle. Referring to fig. 1, we see that e = posina/ cos a, which represents 
the distance of the junction to the central axis, coincides with choice ([10|) . 
For the particular solution (|^) we have a sperical cap of radius e. 

In the generic case, h is characteristic of the form of the cap and e of its 
size. This explains the choice of form ( p^T]) of / which allows us to wedge 
in the cone a cap of a given form of any size. In short, the 2-surface of 
the Euclidean 3-space will be called a smoothed cone of half angle a. The 
corresponding spacetime will be said to have a smooth conical point. We 
emphasize that we can choose Iq and e as little as we want thus making the 
smooth cone as sharp as we want, but keeping the quotient lo/e fixed. 

Finally, we shall also require the continuity of the second derivatives of 
the metric at the junction. This is simply achieved by putting 

r(/o)=0 or h"{^j)=Q. (13) 

As a consequence the Ricci tensor is continuous at the junction. Let us 



note that the supplementary condition (|T^) is not verified by the particular 
solution (y). This was obvious from the beginning since the energy density 
cr is constant in the interior of a straight string and decreases sharply to zero 
at the boundary. The supplementary condition (|TB|) is physically natural 
for a thick relativistic cosmic string but not absolutely necessary. Its real 
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Figure 1: Smoothed cone 
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usefulness is that the continuity of the Ricci tensor makes easier further 
calculations in the neighborhood of the boundary between the cap and the 
cone. 

Ill Self-gravitating string of arbitrary shape 
and coordinate system 

In the previous section we described a straight self-gravitating string with 
some thickness. The aim of this section is to extend this construction to 
a string of arbitrary shape. If the string thickness is sufficiently small and 
if the central line of the string spans a sufficiently smooth world sheet, we 
can suppose that the string approaches a straight one and thus it can be 
locally characterised by a smooth conical point. This suggests to determine 
the spacetime in the neighborhood of a small portion of the string in the 
following way : 

- the interior metric given by 

ds'jj^T = 9AB{r^, I, (j))dT^dT^ + dl^ + f{l)d(j)'^ < I < k (14) 

where /(/) = th{l/t) was specified in Sec. II and the metric components qab 
(y4, 5 = 0,3) are smooth. 

- in the vicinity of the string, the exterior metric given by 

ds%xT = 9AB{r^,l,(t))dT'^dT'^ + df + sm^a{l-lQfd(t)^ I > k- (15) 

We note that a priori we have ommited cross terms in (0) and ([T5|) in order 
to simplify the calculus. Metric (|T^ is regular but the coordinate system 
breaks down at / = 0. We introduce coordinates p"' {a = 1, 2) which are well 
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defined 



= I cos (f) and = I sin (j). 



(16) 



Then, metrics and (|T5|) become 



ds"^ = gAB{,r^.p'')dT 



■^dr^ + gab{p'')dp-dp\ 



(17) 



Let us introduce the function 



r(/) 



with / = y^(pl)2 + (p2)2 



(18) 



/4 Z2 



which is well defined and smooth in the interval < / < Zq from (|^). It has 
null odd derivatives at / = 0. The components gab have the simple expression 



gn = l + r{W)\ gi2 = -r{l)p'p\ ^?22 = 1 + riWf (19) 



in the interval < / < /q- It is interesting to note that gab = Sat + 0{p^). 
There will be no need in the following to express explicitely gab for I > Iq. 

The central line of the string defined by / = spans a timelike worldsheet 
parametrised by whose induced metric is simply 



A radial geodesic on the 2-surface = constant, is also a geodesic of the 
spacetime, normal to the worldsheet at the point P{t^). The coordinate I 
represents the length along this geodesic measured from the point P(r^). So, 
the 2-surface = constant is generated by the spacetime geodesies tangent 
at P{t^) to a 2-plane orthogonal to the worldsheet. Hence we recognise in 
the coordinate system known system (see for example [0) where 

are geodesic coordinates and from which it is easy to extract the extrinsic 
curvature KaAB of the worldsheet by expanding the metric components 



lAB{r'')=gAB{r^,p'' = 0). 



(20) 



gABir^'^p-) = 7AB(r^) + 2KaABir^)p'' + O^p"). 



(21) 
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As we shall see in Sec. V, the extrinsic curvature is implicated in the dy- 
namics of a self-gravitating string. 



IV Expansion of geometrical quantities 

The aim of this section is to expand in powers of 1/e geometrical quantities 
such as the metric, connection and Ricci tensor when e, that is a measure of 
the thickness of the string, is small but not null. Let us point out that all 
length parameters p^, Iq, e are of the same order since lo/e is constant. 

In the interval < I < Iq all these quantities depend on the function r(/) 
through the metric components gabip"")- However since r(/o) = 0(l/e^) it is 
preferable to substitute r(Z) by a function of Z/e 

Qil) = Ml) = '^h'ii) - ^. (22) 

It is evident that q{lo/e) is fixed and depends only of the angle a. The same 
holds for its first and second derivatives q'{h/e) and q"{lo/e). 

The metric components (0), their inverse, and the determinant can now 
be written 

9ab = Sab + q{-)e:e,'^^, (23) 



= > (24) 



A I 

9= detigab) = 1 + ^g(-) (25) 



where e^^ is the totally antisymmetric Levi-Civita symbol and ( ) stands for 
the induced geometrical quantities of the 2-dimensional smoothed cone. 
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Then one immediately gets 

(^a6)^=/o = 0(1), {g'^'U, = 0(1), {9U, = 0(1). (26) 
By direct calculation of the first derivative one obtains 

{9ab,c)l=l, = (g^^^^,^ = 0{-l (27) 

It is also easy to see that 

{9Ab)i=Io = 0(1), {9AB4)l=lo = 0(1), {9AB,c)l=lo = 0(1). (28) 

Now we can calculate the order of the non null connection coefficients for 
l = lo 

Tbc = l9^''i9DB,c + 9dc,b - 9bc,d) = 0(1), (29) 

^aB = \9^''9AB,a = 0(1), (30) 
^B = -\9'''9AB,b = 0(1), (31) 

= ^^'='(^da,6 + 9d,,a - 9at4) =f 1 (p") = 0(1). (32) 

We can now deal with the Ricci tensor for I = Iq. The derivatives of (|3TD will 
give terms of the order of 1/e, whereas the derivatives of (p7[) and (|32D will 
give terms of the order 1/e^. In obvious notations we obtain 

Rab = Rab{-) + Rab{X) (33) 

in which we have 

i?AB(^) = Wb + T^A, (34) 

-Rab(I) =RaB +^aD^"AB " ^^B^'aC ~ ^"BD^^a (35) 
with Rab= do^AB - ^B^AD + ^CD^SiB - ^BD^AC- 
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For the small case indices we get 



Rab = Ral>{\) + Rab{-) + Ral>{l) (36) 



e2' e 
in which we have 

Rab{\) =Rab (37) 

where Rat is the Ricci tensor associated with the connection (^) of the 
smoothed cone, 

Rab{-) = ^dD^ib (38) 

and 

Rabil) = -dtT^, - Tf^T^c- (39) 
Finally for the mixed indices we get 

RaB = -RaB(l) = QaT^j^ - ObT^^ + F^^Ff^ - T^c^aD' (^O) 

The curvature scalar R for I = Iq is given by the equation 

R = R{\) + R{-) + R{l) (41) 



where 



1 . A „^ A 



Rii) =Rab 9^' =R, (42) 



Ri]) = daT%g^^ + T:,T%g^^ + Ti,T^,g^', (43) 

=R +r^D^lB9'''' - rS^r^^^;^^ - r-^^Fl^^^ - d.r^^g'^' - rf^vE^g'^'' 

(44) 

with R=Rab g^^- 

We have calculated the above quantities for / = Iq since we need these 
estimations precisely at the boundary in order to obtain the equations of 
motion in the following section. 
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V Equations of motion of the string 



We now focus attention on the self-gravitating string of arbitrary shape de- 
scribed by metric (^7\ ) in the coordinate system {t^,p°') where the function 
h is specified in Sec. II. We suppose that the exterior spacetime is a vacuum 
solution : Rap = 0, (a = A, a) for / > Iq. 

The energy-momentum tensor Tap of the extended string is the source of 
the Einstein equations 

Ral3 - ^ga/sR = SnGTap , < I < Iq. (45) 

As a consequence of the matching conditions adopted, the Ricci tensor must 
be continuous at the junction I = Iq. So, the interior Einstein equations (^) 
coincide with the vacuum Einstein equations at / = /q 

(Rapho = 0. (46) 

These boundary conditions impose some constraints on the timelike world- 
sheet swept by the central line of the string. We shall examine (|46|) when 
the parameter e is arbitrarily small. 

According to (p6D-(p^), the components (a, b) of equations ( ^6] ) are written 



{Rat)i, = ir^n)ioi^tX + (Ratmo = 0, (47) 

equation being obtained by noting from (^) that [Rab{^)]io = {Rab)io = 
since the junction / = Iq belongs to the cone. 

We need to calculate the limit of the connection coefficient appearing 
in equation (|47|). From (|30|), (pO|), (|21|) and since lo/e is constant we obtain 

limiT^X = lim (Fg,),, = ^^^K,ab = K, (48) 
14 



where Kd is the mean curvature of the timehke worldsheet. On the other 
hand we know from (32) that iXtblio order 1/e. Thus, the hmit e — > 



of equation (^Tf ) muhiphed by e gives 

F;^f,{k^)Kd = with k'' = ^ (49) 

where 

Kikn = lim[e(r^.)/o] (50) 

€— »U 

is finite. Equation (|49|) is independent of e or Iq and depends only on the 
azimuthal angle on the boundary or equivalently as indicated of the unitary 
2- vector k"-. 

In the same way using (|33| ) -(|35| ) along with (|3T| ) ,(|3^ ) , we obtain the zero 



limit of e for the {A, B) components of equations (^) 



where 



F\k'')KbAB = with P = ^ (51) 



F^(r) = lini[e(9,^?^^ + T:,g''\] (52) 



is finite. By introducing the quantities independent of e 



A = g(^)4 and B = q'{^-^A (53) 



after some algebra we can express F (k^) in the following way 



Let us note that the mixed components (a, A) of (^Gf ) give no constraint. 

If the coefficients F^l(fc'') and F^(fc'^) are non vanishing, which is the 
general case, equations (^) and (|51|) yield respectively 

= 0, (55) 
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KiAB = 0. (56) 

Equation ( |5B| ) means that, at the zero hmit of e, the worldsheet swept by 
the central hne of the string is totally geodesic. This of course implies (|55|). 
We recover the situation described in the litterature ^, H when the 
string is a self-gravitating singular line. 

However we see from (|5^) that if we take the condition 

B = -4A, (57) 

then one annihilates the extrinsic curvature coefficient and thus constraint 
( |56D disappears. Hence we are left with constraint ( pS]) which expresses that 
in the zero limit of e the world sheet is minimal. In other words the worldsheet 
has the same evolution Nambu-Goto string. 

We must verify however that when condition (|57|) is applied then equation 
does not also disappear. Equations (EDI) are explicitely expressed as 



e[P + Q{k'f]K^ + k^[P + Q{k'y]K, = 0, 
-k^QiPfK, - e[2P + Q{k''f]K2 = 0, (58) 
-k^[2P + Q{k^f]Ki - k^Q{k^fK2 = 

for all A;^ where P = B + 2A and Q = ~B + AB + \A^ . If A and B 
verify condition (^), then the three equations (|58|) are reduced to only one: 
k^Kx — k^K2 = for all k"'. Hence ( ^5|) is again verified. 

By their definition (|5B|), the quantities A and B are independent of the 
length parameters and depend only on the angle a. That is, we can express 
condition ( ^Tf ) in terms of the angle a by using (|l^) and (|2^) . We obtain the 



simple expression 

sin a cos a = — . (59) 
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This relation is a supplementary constraint on h as it will be explained in 
the conclusion. 



VI Energy-momentum tensor of the string 

The energy-momentum tensor Sap of the string is obtained by an integration 
of Ta/3 on the section = constant of the string. 

Saf^ = f Taf.'/^dpUp^ (60) 

where now e and thus /q are small but fixed. 

From the Einstein equations (|i3| ) and the algebraic expressions of the 
Ricci tensor given in Sec. IV, we can write in obvious notations 

T„;3(r^, p'') = Tap{^) + T,^(^) + T,^(l) / < lo. (61) 

We first apply formula (|60|) for the {A, B) components. Only the first term 
of (^TJ) will yield a non null finite integral when e tends to since the volume 
element of the smoothed cone is of the order /q. 

Sab= f TAB{\)'f9dp'dp' + 0{t). (62) 

Now we have 

TM^) = ~^9AB R= -^9abK (63) 
where K is the Gauss curvature, therefore 

Sab = -^I QabK^Q dp" dp' + 0(e). (64) 
sttG Ji<io 

Since Zq is small we can also approximate the metric qab by development 
(^. We finally obtain 

Sab = -^7AB f K\f^ dpHp" + 0(e) = -p^AB + 0(e) (65) 
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where /x is the hnear mass density. The Gauss-Bonnet formula gives 

/ K^lidp^dp^ = 27r(l - sin a) (66) 

Jl<lo 

and consequently /i is related to a by formula (H). This result is not obvious 



a priori. It originates specifically in the choice of the form of the metric (|1J) 



where / is given by A general discussion of the problem of obtaining 



the energy momentum tensor for a concentrate distribution of matter can be 



found in |]T4| 



The integral of the {a, A) components obviously vanishes as e^. For the 
integral of the (a, b) components the non null finite part vanishes since Rab — 
^QabR = for a 2-dimensional manifold. Thus the integral vanishes as e. 

Discarding terms in e, Sa/s reduces to Sab given by (|65D which is the 
energy-momentum tensor of the Nambu-Goto string. This last result which 
is proved independently of the equations of motion is in accordance with 
them. It can justify a posteriori the definition of a self-gravitating string as 
a smooth cone with metric ( p!4D and (pISl). 

VII Conclusion 

In this paper we have investigated the dynamics of a self-gravitating string 
in the limit where the thickness e becomes negligible. In the generic case 
the worldsheet swept by the string (more precisely by the central line of the 
string) is a totally geodesic surface (the extrinsic curvature of the worldsheet 
is null). This result could have been somehow guessed since the strings 
defined as singular lines of conical points have precisely this property ^ 
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However we have found another pleasant possibihty : by imposing relation 



(59) the extrinsic curvature is no more relevant and only the mean curvature 
is null. This expresses that the worldsheet is extremal which is the behaviour 
of the Nambu-Goto string. If we suppose that the linear mass density fi is 
given, then the angle a is fixed by equation (pBI). It is easily seen that relation 
( |59D is a constraint on the function h; it can be rewritten 

h{^j)h'{^j) = ^. (67) 

The above equation is an algebraic relation taken at the fixed point lo/e which 
is added to the matching conditions (|12|) and (0). It is easily shown that 
such a function h can be found in a polynomial form. The simplest solution 
is an odd polynomial (cf. (4)) of order seven : 

K^) = ^ + &3(^)^ + &5(^)^ + MV- 

see e e 

The four unknown quantities Iq/^, &3, &5 and 67 can be determined as functions 
of the angle a (or equivalently as a function of the linear mass density /i of 
the string) by the four equations (IH), (0) and (^). Of course by taking 
a polynomial of larger order one can get an infinity of solutions. It would 
probably be interesting to investigate the physical meaning of the constraint 
( |67D on the matter of the string. 

By using the method described in Sec. IV, we can evaluate the Riemann 
tensor {Rai3-ys)io at the junction / = Zq or equivalently the Weyl tensor since 



the Ricci tensor (^) vanishes. This gives the magnitude of the gravitational 
field near the thin cosmic string. In the generic case where the extrinsic 
curvature KaBC vanishes, the Riemann tensor is bounded as the radius Iq 
becomes arbitrarily small. On the contrary when only the mean curvature 
Ka is null, the Riemann tensor components {RaBCdjio blow up as 1/e, i.e. 
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l//o, when the radius Iq tends to zero. We see that the gravitational aspect 
is completely different in the latter case. Let us add that if we require that 
the Riemann tensor remains bounded then we necessarily fall in the generic 
case. 



20 



References 

[1] A. Vilenkin and E.P.S. Shellard, Cosmic strings and other topological 
defects (Cambridge University Press, Cambridge, 1994). 

[2] W. Israel, Phys. Rev. D 15, 935 (1977). 

[3] A. Vilenkin, Phys. Rev. D 23, 852 (1981). 

[4] J.A.G. Vickers, Class. Quantum Grav. 4, 1 (1987). 

[5] V.P. Prolov, W. Israel, and W.G. Unruh, Phys. Rev. D 39, 1084 (1989). 

[6] W.G. Unruh, G. Hayward, W. Israel, and D. McManus, Phys. Rev. Lett. 
62, 2897 (1989). 

[7] C.J.S. Clarke, G.F.R. Ellis, and J. A. Vickers, Class. Quantum Grav. 7, 
1 (1990). 

[8] W. Israel, in Journees Relativistes 1989 (Tours) edited by C. Barrabes 
and B. Boisseau, Annales de Physique 14, CoUoque n° 1 (1989). 

[9] L. Marder, Proc. Roy. Soc. London A 252, 45 (1959). 

[10] J.R. Gott III, Astrophys. J. 288, 422 (1985). 

[11] W.A. Hiscock, Phys. Rev. D 31, 3288 (1985). 

[12] B. Linet, Gen. Rel. Grav. 17, 1109 (1985). 

[13] B. Boisseau, and P.S. Letelier, Phys. Rev. D 46, 1721 (1992). 

[14] R. Geroch and J. Traschen, Phys. Rev. D 36, 1017 (1987). 



21 



